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摘要 : 多 孔 介 质 间 突变 界面 在 油气 藏 开发 中 广泛 存在 。 传 统 的 渗流 理论 认为 压力 和 达 西 速度 在 非 连续 界面 处 是 连续 
的 ， 我 们 在 以 往 的 研究 对 其 进行 了 证 仿 。 多 孔 介质 物性 突变 界面 为 固定 界面 ， 本 文 在 前 期 不 可 压缩 流体 的 基础 上 ， 
以 气 - 水 两 相 为 例 ， 考 虑 流体 的 压缩 性 ， 建 立 了 这 类 突变 界面 处 流体 压力 和 渗流 速度 的 连接 条 件 。 研 究 结 果 表 明 , 在 
多 孔 介质 固定 的 物性 突变 界面 : (1) 总 体 压力 必然 连续 ， 但 各 相 流体 压力 可 以 不 连续 (20 总 渗流 速度 相等 ， 但 
各 相 流体 的 渗流 速度 可 以 不 连续 ; (3) 固定 突变 界面 与 移动 突变 界面 的 连接 条 件 不 同 , 不 能 用 Rankine - Hugoniot 
条 件 描述 ， (4) 根据 突变 界面 的 连接 条 件 以 及 两 相 达 西 公式 根据 界面 两 端 压力 以 及 流体 分 布 可 计算 出 跨 物 性 突变 
界面 的 各 相 流 体 的 流速 ， 跨 物性 的 突变 界面 的 流量 应 取 上 游 值 。 
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Jump interface in porous media exists widely in oil and gas reservoirs during development. The traditional seepage 

theory holds that the pressure and Darcy velocity are continuous at the discontinuous interface, and we have 

falsified them in previous studies. In this paper, taking gas-water phase as an example and considering the 
T= compressibility of the fluids, the connection conditions of phase pressures and Darcy velocities at the immovable 
> jump interface are established. The results show that :(1) the total pressure is always continuous, but the phase 
©) pressure of each phase can be discontinuous;(2) The total Darcy velocity at the both sides are equal, but the 
©9 phase Darcy velocity probably discontinuous;(3) The connection conditions of fixed jump interface do not comfort 
= to Rankine - Hugoniot conditions, which are different from mobile jump interface, i.e. flood front which 
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= 突变 界面 在 油气 藏 中 广泛 存在 ， 例 如 储 层 的 层 间 差异 、 原 始 的 油气 水 界面 、 生 产 过 程 中 注水 或 注 气 前 缘 、 化 

L 学 驱 油 前 缘 、 边 底 水 锥 进 界面 、 基 岩 与 裂 颖 系统 的 国 连接 界面 、 井 简 与 储 层 的 连接 面 、 压 裂 颖 与 储 层 的 交界 面 ， 

O 相关 的 多 相 渗流 问题 在 数学 上 称 为 Neumannn 问题 。 两 相 渗流 突变 界面 处 的 连接 条 件 是 渗流 数学 模型 的 重要 组 
成 ， 是 时 空 行 流体 界面 跟踪 ( 示 踪 剂 )、 计 算 油 气 并 产量 、 描 述 跨 类 型 或 跨 尺 度 多 和 孔 介质 流体 渗流 规律 、 相 渗 测 试 、 
构建 油 藏 数值 模型 等 多 相 渗流 领域 的 基础 理论 问题 。 根 据 渗流 文献 ， 在 多 孔 介质 突变 界面 处 流体 各 相 压力 和 渗流 
速度 连续 。 即 : 


Vy = XV =V; 
D, = Py Dy B, (1) 
p, =P, 
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根据 质量 守恒 方程 ， 单 位 时 间 内 微 元 内 流体 变化 F 流 入 微 元 的 流量 流出 微 元 的 流量 。 对 于 移动 界面 ， 体 积 关 
0， 没 有 源 汇 项 ， 则 有 ， 了 流入 微 元 的 流量 | 流出 微 元 的 流量 F0。Hassanizadeh, S. M., & Gray, W. G. (1989) 用 物质 守恒 
方程 对 其 进行 了 证 明 。 针 对 界面 两 侧 区 域 毛 管 力 不 同 的 情况 ， 作 者 在 文献 中 对 传统 的 界面 条 件 进行 质疑 ， 然 后 根 
据 物 质 守 恒 方程 和 两 相 渗 流 公式 推导 了 新 的 界面 条 件 以 及 跨 突变 界面 的 渗流 公式 路 习 。 气 - 液 两 相 驱 蔡 界 面具 有 三 个 
特征 : COD 多 孔 介 质 物性 突变 界面 不 随时 间 不 断 移动 (2) 在 多 孔 介 质 的 突变 界面 处 流体 分 布 可 能 不 连续 ; 
G) 气体 压缩 性 高 ， 驱 蔡 过 程 需 要 考虑 压缩 性 。 所 以 本 文 将 在 前 期 不 可 压缩 流体 两 相 突 变 界面 条 件 的 研究 基础 
上 ， 考 虑 流体 的 压缩 性 ， 以 气 - 水 两 相 为 例 ， 研 究 多 孔 介 质 突变 界面 处 ， 气 - 液 两 相 压 力 和 渗流 速度 的 连接 条 件 。 


2 传统 突变 宜 面 条 件 面 对 的 挑战 中 
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CL). 突变 界面 的 连 扎 条 件 针 对 界面 两 侧 不同 的 质点 ， 建 立 的 是 不 同 质点 压力 的 连接 条 件 和 不 同 滩 流 速度 的 连 
接 条 件 〈 如 图 D ， 该 条 件 所 描述 的 并 非 突变 界面 上 同一 流体 质点 的 流速 和 压力 。 传 统 的 渗流 理论 认为 界面 体积 
0， 根 据 质量 守恒 方程 得 出 脏 入 突变 异 面 的 流量 等 于 流出 突变 春 面 的 流 员 这 一 结论 虽然 是 正确 的 ， 但 它 描述 的 并 全 
界面 的 条 件 。 如 图 1 所 示 ， 以 一 维 为 例 ， 假 如 包括 界面 的 薄 层 状 微 元 厚度 为 六 9， 流 入微 元 的 流速 记 为 v7 ， 流 出 
微 元 的 流速 记 为 v* (4 = 0,w) ， 两 者 的 数量 关系 即 为 渗流 速 的 连接 条 件 ， 应 满足 根据 质量 守恒 方程 2)， 根 据 该 广 
程 并 不 能 得 出 V; = vi: 而 传统 的 渗流 力学 针对 界面 条 件 所 用 的 质量 守恒 方程 表达 为 (3)。Hassanizadeh, S. M., & 
Gray 试图 证 明 由 (2) 可 以 得 到 (3) ， 但 我 们 发 现 证 明 过 程 有 误 。 
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图 1 界面 条 件 与 界面 上 质点 守恒 方程 所 指 对 象 图 2 2 包括 驱 痊 界 商 的 江洲 层 状 微 元 示意 图 


(2) 研究 一 维 水 驱 油 的 Buckley-Leverett 方程 的 解析 解 显示 渗流 速度 在 前 缘 界 面 处 不 连续 分 布 ， 这 与 传统 的 
界面 条 件 〈 如 公式 (1D)) 明显 矛盾 。 此 外 De Neef, M. J., & Molenaar, J. (1997) 针 对 突变 界面 两 侧 毛 管 力 曲 线 不 相同 ， 
对 传统 的 界面 条 件 提出 了 改进 ， 他 们 认为 润 相 流体 压力 连续 ， 但 非 润 湿 相 压力 可 以 不 连续 ; 

G) 调研 文件 我 们 发 现 虽然 (1) 表 示 的 界面 条 件 广泛 见 诸 文 献 ， 但 并 没有 得 到 真正 应 用 ， 因 为 公式 (1) 所 示 的 传 
2 件 ， 以 及 两 相 流 达 西 公式 ， 突 变 界 面 处 传导 率 应 该 为 界面 两 侧 的 调和 平均 值 ， 它 不 能 正确 反映 真实 的 
参 流 过 程 ， 所 以 被 单 点 上 游 权 (SPU)〉 和 其 它 高 阶 上 游 权 方法 取代 。 
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(4) 从 计算 角度 ， 即 使 在 突变 界面 处 ， 流 体 分 布 昌 然 突 变 但 仍然 连续 的 情况 下 ， 由 于 突变 区 分 布 厚度 薄 ， 将 
突变 区 简化 为 不 连续 界面 ， 可 实现 在 获得 相近 计算 结果 的 前 提 下 ， 大 幅度 减 小 计算 量 。 


0 T : 
图 3 GREN EGG T"— 卖 界面 示意 图 
CE: 真实 界面 附近 流体 分 布 ， 右 :简化 后 的 不 连续 界面 ) 


3. 考 虑 突变 界面 的 两 相 可 压缩 流体 Darcy 公式 


3.1 质量 守恒 方程 

假设 油气 藏 存在 突变 界面 ， 界 面 位 置 不 随时 间 变 化 ， 界 面 两 侧 储 层 的 孔隙 度 ， 渗 透 率 ， 流 体 压力 ， 流 体 饱和 
度 ， 孔 隙 度 ， 着 石 压缩 性 不 同 。 整 个 油 藏 被 突变 界面 分割 为 两 个 子 区 域 ，@O AO, (图 2)。 每 个 子 区 域 质量 守恒 
方程 满足 微分 方程 : 


o 
Fls) = Sy (pva) 
(5) 
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式 中 下 标 《 为 流体 类 型 标记 ， 可 选 气相 和 水 相 ; ØRSTA: p TAC NERE Kg/m*; V 表 示 渗 流速 度 ， 
mid; t 表 示 时 间 ,d。5S, 表示 流体 饱和 度 。 两 种 流体 的 饱和 度 有 如 下 关系 : 


5,+58,=1 (6) 
引入 地 层 体积 系数 的 定义 公式 : 


pote p -2e (7) 


其 中 已 表示 气体 在 地 层 条 件 下 的 体积 与 地 面 气体 休 积 的 比值 ， 亦 可 表示 气体 在 地 面 的 密度 pp CHY kg/m) 与 地 
层 密度 o, (单位 kg/m) 的 比值 ，B, 表示 水 在 地 层 条 件 下 的 体积 与 地 面条 件 也 的 体积 ， 亦 可 表示 水 在 地 面 的 密度 
pus HUE kg/m) SEFE p, CHY kg/m) 的 比值 。 方 程 组 C5) 改写 为: 
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将 上 述 方程 组 左 端 关于 时 间 的 导数 项 展开 ， 并 引入 岩石 压缩 系数 cy 〈 单 位 MPa’) ， 气 体 压缩 系数 c。( 单 位 MPa’) 


O 以 及 水 压缩 系数 C， CBE MPa O ， 方 程 组 改写 为 ; 
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3.2 储 层 物性 突变 界面 处 的 渗流 速度 的 连接 条 件 
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将 程 组 〈9) 中 气相 方程 和 水 相 方程 分 别 除 以 饱和 度 导 数 项 的 系数 ， 然 后 相 加 得 到 两 相 整 体 流动 方程 


或 


ap op d (v d ( v 
S, (c; *6,]—7 9S, (c; +c, ]J— - -B 3. |-B - (10) 
$5. (c, 2 ôt #8, (c, ;) ôt dix Ae] 
P 气 、 水 两 相 流 体 的 渗流 满足 两 相 达 西 公式 : 
vi o7 KA [S5 pgsina) (4^2 w,g) AD 
dx 
OP nodu ve (=w,g) (12) 
dx KA, 


T" k,, E 
4 表示 流 度 ，4 ==, ， 上 ,为 流体 《的 相对 渗透 率 ，j RANT C BUR, mPas; pA CR, MPa: 
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对 方程 (10) 两 侧 在 [x , x 区 域 关 于 x 做 积分 有 


Xp Q Xp 
tete te eee p B, AE a (se (13) 
g w 


[xo x,] 区域 包括 突变 界面 ， 所 以 方程 (13) Æ GUN LHS) 有 ; 


a o 
LHS - (x, se «Tees; (c; m: 


(14) 
e b 
+(x, ne [sug vet) sse) 
MAINA : 
d| v dp . V 
By |=v +v c sy +vc sin a -— (15) 
dx g 99 dx g [as & 
d [v dp : V 
B —|  /=v.4+vic “=y +V C sin = (16) 


将 (15) (18) 代入 方程 (13) 右 端 〈 记 为 RHS) 各 项 ， 于 是 : 


pe 20 ¥) ys, 门 - ina- ing - s an 
X SIN C 一 SVC sin X 
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上 式 积分 项 内 所 有 参数 都 是 真实 世界 中 数量 有 限 的 物理 量 ， 所 以 有 积分 后 有 : 
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由 (4) 和 (18) 得 到 如 下 方程 : 


(x, —x,)9* se «cas +S° ‘(os x 


(18) 
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cse siete) 7 +S2 (c ec) — 
" (19) 
=vi va +(x, —x,)| vc | påg sina — js lap påg sina — ur 
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又 因为 x —xo xx MWA 


vf =y” (20) 
上 式 说 明 ， 在 固定 的 突变 界面 处 ， 气 和 水 的 总 的 渗流 速度 的 分 布 连 续 。 
3.3 储 层 物性 突变 界面 处 的 压力 的 连接 条 件 
由 气 、 水 两 相 的 达 西 公式 〈11) 得 到 : 


dp dp. . 
,=w+V =-K4, ES f. 2. (f,0,9 + f. o. g)sina (21) 


定义 总 体 压力 (total pressure? 


1 . Å 
D=p,+f, f. de dS, (22) 


1. dp. 1. dp. 
dp, do of fu Sw om 6j. f. ge as, 
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dx dx oo Ox OS Ox 
定义 办， 因为 气体 粘度 随 压 力 变化 平缓 ， 工 程 上 可 采用 如 下 公式 近似 
1 dp 1 dp 
O € dS o < dS 
f fy dx w 7 N 人 dx W 


OD D Op, 
将 (24) RA (22) 可 以 得 到 : 
d 
Li 09 e oig abs 
dx dx dx 


将 其 代入 〈21) 总 体 流速 表达 式 得 到 : 


整理 得 到 : 


两 端 分 别 对 x 做 积分 有 : 


p2a- er rre a D 

% dx « -K4 (1 
积分 后 可 得 
oo- | AA , No | ee (6-7 
o7 KUAMI-o") K'a*(1- o") 1- o* 1-0" 
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上 式 说 明 在 突变 界面 处 有 总 体 压 力 连续 。 

3.4 物性 突变 界面 处 的 流体 压力 和 渗流 速度 的 连接 条 件 

根据 〈20) GO 可 得 到 多 孔 介 质 在 固定 的 突变 界面 处 流体 压力 和 渗流 ; 
fi (v: -GPC*) = fs (vs -GPC") 


f; (vi -GPC")= få (vi +GPC") 


1 dp* 1 dp? 
M ii pf. pe (ey PS) ae 


dg 
ng 4j dr G) 
e eE) ag p- clo} 
p;-[.. fz( E=pr-få F RE 
A; AT a ALAP b 
其 中 : GP! = 天 "一 一 | Api gsi , GP’ = 天 "一 An 
AL dx T dx 


根据 (310 可 知 ， 在 物性 突变 界面 处 流体 压力 和 渗流 速度 可 以 不 连续 分 布 。 
4 跨 突变 表面 的 流量 计算 


根据 (29) 可 得 : 
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速度 的 连接 条 件 : 


(31) 
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中 pat DL x)? (x, x)? 
v = l-@ 1-o (32) 
Xp — X, X, — X. 
KX (1-w") ou o") 
将 (320 RA (QD "EZ HSE vi v; 
a a a dp: dp: (S;) s) OS? a a : 
Vw fo t ds? Ox +(ps—ps)gsina 
(33) 
V fv dp: ( Jos; ( å sina 
a as ae ME 
LAVE, v 
VENE PUSSE å ak. 
v= fv + ds?! ex «(os p; gsina 
P (80) : (34) 
Pe (9w) OS? . 
ve = fev’ as) ox (e p,)g sina 
2 跨 物 性 的 突变 界面 的 流量 应 取 上 游 值 ， 即 : @ 当 流体 由 Xx, 流向 XxX,， 通 过 突变 界面 的 流速 为 vs vo: QUU 
2 ` HH x, DID x,» JEGERE FET AEE v^ ， vi 
REDI. 
| (OD 多 孔 介 质 物 性 突变 界面 是 一 种 固定 界面 ， 界 面 处 流体 压力 和 渗流 速度 的 连接 条 件 可 表述 为 : 总 体 流速 相 
等 ， 总 流速 相等 ， 但 流体 压力 和 渗流 速度 可 以 不 连续 。 流 体 滩 流速 度 不 连续 是 因为 渗流 速度 记录 的 是 流量 除 以 渗流 
Se 模 截 面积 ， 即 使 流体 分 子 的 真实 速度 在 界面 处 连续 ， 但 渗流 速度 可 以 不 同 ; 微观 上 当 流 体 在 界面 处 连续 分 布 ， 流 体 
O 压力 场 在 界面 处 应 连续 , 但 在 突变 界面 所 处 , 两 相 流体 很 容易 在 孔 喉 位 置 发 生 截 断 , 从 而 形成 压力 场 的 不 连续 分 布 ; 
人 (2) 对 于 突变 界面 处 流体 是 否 能 够 不 连续 分 布 ， 学 术 界 还 有 争议 ， 一 些 学 者 根据 传统 的 界面 条 件 认 为 流体 不 
E 可 能 连续 分 布 。 我 们 认为 ， 抛 开 争 议 ， 从 等 效 计算 的 角度 ， 流 体 即使 在 突变 区 连续 分 布 ， 但 分 布 范围 很 狭窄， 渗流 
(5 阻力 占 比 极 低 ， 采 用 本 文 提出 的 能 够 反映 跨 间 断 界 面 渗流 规律 的 界面 连接 条 件 ， 可 避免 采用 高 密 网 格 采 样 条 件 下 ， 
得 到 满足 计算 结果 ; 
GO 大 多 数 情况 下 ， 分 相 流 系数 随 压 力 变化 较 小 ， 可 忽略 C@=0) ; 
(4) 多 孔 介 质 物性 突变 产生 的 界面 peru ree 与 相 比较 ， 流 体 渗流 速度 和 压力 在 界面 的 连接 条 件 
特征 不 同 ， 即 静止 的 突变 界面 的 连接 条 件 与 移动 的 突变 界面 的 连接 条 件 不 同 ， 需 要 注意 应 用 条 件 。 例 如， 本 文 所 建 


立 的 界面 条 件 不 适合 描述 跨 尺 度 ， 跨 类 型 介质 的 流量 的 计算 ; 
(5) 本 文 以 气 -水 两 相 为 例 ， 没 有 考虑 溶解 性 ， 如 果 考 虑 溶解 性 ， 
之 变化 ， 本 文 只 考虑 了 一 维 情况 ， 对 于 多 维 界 面条 件 更 为 复杂 。 
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